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Abstract. The eigenvalues of the eUiptic A'^-body Ruijsenaars operator are obtained by a dynamical 
version of the algebraic nested Bethe ansatz method. The result is derived by using the construction 
given in |Q, where the Ruijsenaars operator was obtained as the transfer matrix associated to the 
symmetric power of the vector representation of the elliptic quantum group i?T,7{gljv)- 



1. Introduction 

The goal of this paper is to solve the elliptic A^-body Ruijsenaars model by the Bethe ansatz method. 
It was done in ||^ for = 2, as a special case of the diagonalization of the transfer matrix of particular 
highest weight modules associated to the elliptic quantum group £^T-,^(gi2)- The result was achieved by 
a dynamical generalization of the algebraic Bethe ansatz (see ^, |^ for the description of the algebraic 
formulation of the Bethe ansatz). For general N , the Ruijsenaars operator can be constructed as the 
transfer matrix T{z) associated to the symmetric power of the vector representation of £'r,7(gljv) [Ql- 
Since Tiz) is a trace over an iV-dimensional space, we diagonalize it by the nested version of the Bethe 
ansatz 

The Ruijsenaars operator is a difference operator which is a q-deformation of the Calogero differen- 
tial operator. In the solution of the elliptic A^-body Calogero model (in the Bethe ansatz form) was 
obtained from the integral representation of solutions of the elliptic Knizhnik-Zamolodchikov-Bernard 
equations by applying the stationary phase method |@, § . Similarly, there is a link between the Rui- 
jsenaars model and the qKZB difference equations. The eigenfunctions of the trigonometric Ruijsenaars 
operator are given in |0 in the form of Etingof-Kirillov traces of intertwining operators . As for 
the integral representation of solutions of the qKZB equations, it is known for SI2 only JT^ . 

The structure of the paper is as follows. In section^, we review the construction of commuting transfer 
matrices associated to the representations of the elliptic quantum group Er.-^(g^i^)- In section ^, we 
recall how the Ruijsenaars operator is related to the transfer matrix T(z) associated to the symmetric 
power of the vector representation of Et^^{^j^). In section we explain the idea of the dynamical 
version of the algebraic nested Bethe ansatz. Finally in sections || and we apply explicitly the Bethe 
ansatz method; we write down the Bethe equations and the eigenvalues of the transfer matrix T{z). 

2. The elliptic quantum group associated to gl^ 

We recall the definition of the eUiptic quantum group Er,'y{g\iq), following Let g be a Lie algebra, 
f) a Cartan subalgebra of g and [)* its dual space. Let 11^ be a finite dimensional diagonahzable f)-module, 
i.e. a complex finite dimensional vector space with a weight decomposition W — ©^g(,*VF[/i] such that 
f) acts on W[ij] by xw = ii(x)w for a; 6 f) , w £ 
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The starting point to define representations of elliptic quantum groups is an i?-matrix R{z, A) G 
End(VK (8> W) depending on two parameters z G C, A G f)* and solution of the dynamical Yang-Baxter 
equation 

(1) -Z2,X~ ih^'^) - Z3, A) i?(23)(^^ _ ^ _ ^^(D) 
= _ ;,3, A) R(''Hz, - Z3, A - jh^'^) R^^'Hz, - Z2, A). 

In this equation 7 is a generic complex parameter, and we use the standard notation 

(2) /(A - 7/1) w = /(A - 7yu) w if w G VF[/Li] 

for any complex- valued function / of A. 

In what follows we are interested in the case g = gl^, and f) the algebra of diagonal complex N x N 
matrix, acting on the vector representation V = of gl^ with standard basis {ej)j^i^,,, ^n- We identify 
I) and ()* via the trace, and with C''^ via the orthonormal basis {luj — ejj)j=i,,,, ^n, denoting ejk the 
N X N matrix acting on the standard basis by Cjkei = Skiej. The weight spaces of the f)-module V are 
V[ujj] — Ccj. In this case the notation (H) reads explicitly 

/(A -7/1) efe = rfe/(A) efe 

where Tk is the shift operator by —jujk'- if / is a complex- valued function of A, 

rfe/(A) = /(A-7cjfe) = /(Ai,... ,Afe_i,Afe -7,Afc+i,... ,AAr). 

Let A) G End(C^ ® C^) be the elliptic solution of the dynamical Yang-Baxter equation 
given by the formula 

JV JV N 

(3) R(z, A) = ^ en ® -I- ^ a{z, \i - Xj) en (g) ejj + ^ Ai - Aj) eij ® eji, 
where 



i=l i,i = l i.3 = l 

i^j 



0{z)9{X + j) eiz + X)9ij) 

9{z--f)e{X) e{z--f)0{x) 

is the Jacobi's first theta function: 



= -^exp 



Mj + ^fr + 27:i{j + ^)iz+^) 



where t is a complex parameter with Im(T) > 0. 9 is analytic with zeroes (n -I- Tm),n,m G Z. It 
satisfies 

9{-z) = - 9{z) = 9{z + l) and 9{z + t) = - 6l(z) e-^'^^+r)^ 
The i?-matrix has the following properties: 

(4) i?(i2)(z,A) i?(2i)(_^^A) = Id, 

N 

(5) i?(i2)(0,A) = P(i2) ^ ^ e,,®e,„ 

(6) A),a;(8ld-f Id® x] = 0, Vx G t), 

(7) [r(i)r(2),i?(i2)(^,A)] = 0, 

where F is the diagonal N x N matrix diag(ri)i=i.... ^n. 
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A representation of the eUiptic quantum group E = Er.-yiglpj) is by definition a pair {W,L) where 
W is a, finite-dimensional diagonalizable t)-modulc and L{z, X) is a meromorphic function with values 
in End(, (C^ (8) W) , obeying the relation 

and commuting with the action of t): 

[L{z,X),x(E)ld + ld'»x] = 0, Vx £ f). 

An i?-submodule of an i?-module {W,L) is a pair {W',L') where W is an ()-submodulc of W such 
that ® W is invariant under the action of L{z, A), and L' is the restriction of L to this invariant 
subspace. £^-submodules are _E-modules. 

The basic example of an i?-module is {C^,L) with L{z,X) — R{z — w,X). It is called the vector 
representation with evaluation point w and is denoted by V(w). 

Other modules can be obtained by taking tensor products: if (Wi, Li) and (W2, L2) are £'-modules, 
then also {Wi (8)^2, L), with an f)-module structure x{wi '8)102) = xwi 'SyiW2 +wi (!S)XW2 and a L-operator 
L(z,A) =L^''^(z,A-7/i(3)) L^^^^^z.X). 

It is useful for what follows to introduce what we call the Lax operator associated to the i?-module 
(W, L). It is an X matrix, with elements which are operators acting on the space of meromorphic 
functions of A G f)* with values in W. It is defined by the formula 

(8) C{z)w{X) = L'^^^Hz,X)T<^^'> w^^\X). 

More explicitly, let us introduce matrix elements by L{z, A) Cj ®w = ^ ^ ® Lij (z. A) w (and the same 
for C). The elements Cij{z) of C{z) act the following way: 

Cij{z) w{X) = Lij{z,X) 'w{X - -luij). 

Using the property of the i?-matrix, we see easily that £ satisfies the commutation relation 

(9) R^'^\z,^Z2,X~^h^'^) C^'^\z,) £(23) (^2) ^£(23) ^(13) R(^^){z,-Z2,X). 

The trace of the Lax operator tr(^^£(-^'^) (z) = Cii{z) leaves the zero weight subspace W[Q] of W 
invariant. The transfer matrix associated to the E'-module [W^L) is, by definition, 



W[0] 



As a consequence of relation (|^), the transfer matrices commute for different values of the spectral 
parameters. 



3. The Ruijsenaars operator 

The elliptic Ruijsenaars operator is (up to a conjugation by a function) a difference operator acting 
on functions of A € I)*: 

M = ^ J]^ fl'/7 tl^^ r,, with a coupling constant £ e N. 

■'=1 

It can be obtained as the transfer matrix associated to a particular i?-module that we now introduce. 
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Let n e N. Let F®"(0) denote the ^'-module V{0) ^(7) V{-f{n - 1)). F®"(0) is the pair 

{W,L) with = (C^)®" = • • • C^, and i is the operator on (C^)®" given by the 

formula 

n n 
3=2 j=3 

We denote by S'"(C^) the space of symmetric tensors of (C^)®", ie. the subspace of (C^)®" 
invariant under the action of the symmetric group 5„. It is proved in Q that S'"(C^) is an i?-submodule 
of y®"(0). It is called the nth symmetric power of the vector representation (with evaluation point 0) 
and is denoted by ^"^(O). The zero weight subspace of this module is trivial unless n is a multiple of 
N; if we take n — N£, the zero weight subspace is one dimensional and is spanned by the sum of the 
tensors e^j (8) • • • SJei^j over all sequences (ij) such that each integer between 1 and N occurs precisely i 
timesFl Let T{z) be the transfer matrix associated to S'^^V{0). If we identify the zero weight subspace 
of ^^^(O) with C, we have 

Thanks to this result, we can solve the Ruijsenaars model by applying the Bethe ansatz method to the 
transfer matrix T(z). 

Let us remark (although we do not use this result in what follows) that M belongs to a family of N 
commuting difference operators 

(where |/| denotes the cardinality of a subset / of {1, . . . , A^}). Each operator M„ can be constructed 
by considering transfer matrices associated to _E-modules obtained as symmetric and exterior powers 
of the vector representation of ET^^{glj^). 



4. Setting up of the Bethe ansatz 

Let us recall that the algebraic Bethe ansatz is a method for diagonalizing transfer matrices obtained 
as a trace of a 2 x 2 Lax matrix with a spectral parameter. The nested Bethe ansatz is a generalization 
of the Bethe ansatz to the case of Lax matrices with an iV-dimensional auxiliary space, for iV > 2; it is 
achieved in (N-1) steps, each step resulting in reducing by 1 the dimension of the auxiliary space. 

Here we apply a dynamical version of the nested Bethe ansatz: the Lax operator is a matrix on 
an auxiliary N dimensional space which elements are difference operators acting on functions of N 
parameters (Ai, ... , Xn) with values in some vector space. 

The idea consists in looking for eigenstates obtained by applying "creation operators" B{tk) to 
a reference state v: ■(/'(^ii ■ • ■ j^m) — B{ti) ■ ■ ■ B{t,n)v. The B{tk)'s are particular elements of the Lax 
matrix C, and the reference state v is chosen to be an obvious eigenstate of the elements of the Lax 
matrix. The problem amounts to finding the conditions on the parameters tk- More precisely, in the 
case N = 2, B{t) is the operator £12 (t). In the case TV > 2, it is a little more complicated: the operators 
B{t) belong to the set {Bj{t) = Cij)j=2.... ,n, so that ijj is roughly of the form Bj-^ {ti) ■ ■ ■ Bj^ {tm)v and 
we have to choose not only the right parameters tk, but also the indices jk- 



^By extension we call zero weight a weight which is a multiple of a; = '}2iLi^i- Since the i?-matrix and the Lax 
operator depend only on differences (A; — Xj), adding the same constant to each A; does not change our results. 
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Let us now go into details. We want to diagonalize the transfer matrix 

(11) T{z) = tr(°)£(2) 
where 

Nt 

(12) C{z) = i?(°i)(z,A-7^/i(^)) ••• i?(°'^^)(z-7(iV€-l),A) T^. 

This Lax operator is a matrix on an auxiliary A'' dimensional space (denoted by the index (0)) with 
elements which are difference operators acting on functions of A = (Ai, ... , Ajv) with values in = 
(C^)®^^. We denote these elements by: 



^^^^ - ^ Ci{z) D,j{z) 



i,j=2,... ,N 

We denote by B the lign vector {B2, . . . , Bff) and by D the matrix (-Dij)i,j=2,... ,jv- 

4.1. The form of the Bethe ansatz. Now let us explain more precisely in which form we are looking 
for the eigenstates of T{z). 

First we choose a reference state v. We take a joint eigenstate of the operators A and Da which is 
furthermore annihilated by all the Dij for i ^ j: 

V = ei (8) • • • (g) ei gW. 
^ ' 

Nt factors 

Proposition. The action of the operators A and Dij on the reference state is given by the following 
formulae: 

(13) A{z) [giX) v] = g{X - 7^1) v, 

N N 

(14) D{z) [g{X) v] = Yl v] = Y^eu Mz,X) g{X-^uji) v, 

i,j=2 i=2 

Where ^ (z A) = ^(^) ejX, - X, + Nlj) 
wnere<s,[z,A) _ ^^^^ ^(A^ - Ai) ' 

Proof. We can write R{z, A) = X]fj=i ^ A) with 

Rii{z, A) = Cji + ^ a{z, Aj - A^) ejj for i = 1, . . . ,N, 

Rij {z, A) = (3{z, Xi - Xj) Cji fov i,j = 1,... ,N, ijt j. 

Then we have 

N Nl Ni 

L,,{z,X) = 4K-,A-7E^''04'U^-^'^-^^'^''' 

ji,---,jNe-i fc=2 fc=3 



Let us consider Ly , for i, j > 2, i ^ j. In each term of the sum we have at least one of the factors 
which is of the form Rik for some k ^ i. Since Rik is proportional to eki, Rik^i = 0. Thus Dij v = 

for i,j > 2, i ^j. 
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Now if we apply La to v, the only term which gives a non zero contribution is the one corresponding 
to ji = • • • = jNt-i = i- Therefore 

A) V = V and Lii{z,X) v = <f>i(z,A) v if i ^ 1, 

Ni 

with$i(0,A) J|a(z-7(fc-l),Ai-Ai+7(7V^-fc)). □ 
fc=i 

Then we look for eigenstates in the form 

N 

ip{ti,... ,t^,X)^ ^ Bj^{ti)---Bj^{t^) V gj,,,,,j^{\), 

with some coefficients Qj^^...^^ which are functions of A (and depend implicitely on the parameters t^). 
We impose conditions on these j„'s so as to consider only the states i/; of weight zero, v has a 
weight [{N — 1)£, —t, . . . , —£] in the basis {luj — Cjj — uj/N). Since applying the operator Bj to a vector 
of weight /i gives a vector of weight ^ — uji + , we have to apply to v exactly i times each of the Bj 
in order to get a zero weight vector. This means that we must take a state of the form 

(15) V(ti,...,V-i)^:A) = i?W(ii)---S«^^^)^Ht(Ar_i),)i;5(i--(^-iW(A) 

where gt^'- - -(^-i)^) is a function of A with values in the zero weight subspace of = (^£^n -lysiN ^ 

4.2. The commutation relations. To evaluate the action of the transfer matrix 

T{z) = A{z)+tTD{z) 

on the vector i/;, we have to push A and D to the right of the B operators, and we need to write the 
commutation relations (^) more explicitly. 

Proposition. The commutation relations can be written in the form 

B^^\z,) B^^\z2) = B^'\z2) B^^Hz,) R^^^Hz,~Z2,X), 

A{zi) B{z2) = B{z2) A{zi) a{z2 - zi, A) + B{zi) A{z2) b{z2 - zi, A), 

(16) D^^^ (zi) ^(2) (z2) = a<-'^ (zi - Z2, A - 7/1) B^^) (^2) D<-'^ (zi) R^^^^ (zi - Z2, A) 

+ cW(zi-Z2,A-7M B^'\z,) i?W(z2) P(i2), 
where A = (A2, . . . , \n), R is the R-matrix associated to gljy^i 

N N N 

R{z, A) = ^ (g) eu + ^ a{z, Xi - Xj) eu (g) ejj + ^ /3(z, Ai - Xj) eij (g) Cj^, 

i = 2 >.3 = 2 i.j = 2 

p(i2) — 'Yl!ij=2 ® ^Ji permutation operator, and a(z, A), 6(z, A), c(z, A) are diagonal {N — 1) x 

(N — 1) matrices with elements {j ~ 2, . . . , N) 



1 1 lu \M Piz,Xi-Xj) /3(z,Aj-Ai) 

Q!(z,Aj — Aij Q!(z,Aj— Aij Q!(z, Aj—Aij 



Proof . If we denote by Rikjn the elements of the i?-matrix: 

AT 

i?'^^^)(z,A) = ^ Rik,]n{z,X) eik'gejn, 
i,j ,k,7i—l 
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the element (ifc, jn) of the matrix relation (|^) reads 

N N 
r,s—l r,s—l 

Taking successively i — j — 1 and k,n > 2; i — j — n — 1 and k > 2; j — 1 and i,k,n > 2, and using 
the particular form of the _R-matrix {Rikjn ~ unless (fc, n) = (i,j)), we get 

JV 

Bk{z) Bn{w) = ^ Bs{w) Br{z) Rrk,sn{z ~ W,X), 
r,s=2 

Bkiz) A{w) ^ Bkiw) A{z) P{z-w,Xi- Xk) + A{w) Bkiz) a{z - w, Xk ~ Xi), 

(3{z -w,{X- 7/1)4 - (A - 7/1)1) Bk{z) Din{w) + a{z -w,{X- -fh)i - (A - 7/1)1) Afe(^) Sn(w) 

N 

= ^ BsC^) Ar(2:) -Rrfc,sn(2 - W, A), 



r,s=2 



which can be written in the form (16). □ 

Since the coefficients appearing in the commutation relations (^6|) depend on A and /i, we also need 
to know the way the functions /(A, h) go through the different operators. 

Proposition. For any complex- valued function f of X and h, we have: 

A{z) f{X,h) = f{X~ju;i,h) A{z), 

(17) B,{z) f{X,h) = f{X--iLo,,h + uJi-oj^) B,{z), for 2,... ,7V, 

Dij{z) f{X,h) = f{X--fujj,h + uJi-ujj) Dij{z), for i, j = 2, . . . , A^. 

Proof . For the A dependence it is a straightforward consequence of the definition (^) of the Lax operator. 
For the h dependence, it is not difficult to see that since L commutes with the action of f), we also have 
+/j(2))j ^ 0, or equivalently C,j{z)f{h + LOj) = /(/i + w,)£„(z). □ 

5. The first step of the Bethe ansatz 

With the help of the commutation relations (16, |l^), the action of T = A + trD on the vector tp 
given by ( |l5|) can be recast in the form 

T{z)^ih,... ,t„,,X) = i3(^'(ti)---i?^"^(t™) z;g^''-'"'(A) 

m 

+ J2B^'\z)B('^+'\tk+^) ■ ■ ■ B^^^\t„,)BW(^t,) ■ ■ ■ B^^^-^\tk-,) V g^^- '^\x). 

k=l 

Here and until the end of this section, we write m = {N — 1)£. Admitting that these different terms 
are linearly independent, we have to impose that gk = for k = 1, . . . ,to, and that go is proportional 
to g. The term associated to go is usually called the "wanted term" and the other ones the "unwanted 
terms" . 

If we carried out the commutation procedure, we would get 2™ terms for the A part, and as many 
for the D part: we see that this is hopeless to apply the process literally for all the terms. Nevertheless 
the "wanted term" and the first "unwanted term" are easy to obtain. 

To get the "wanted term" , we need to keep only the first part of the commutation relations A{z)B{t) = 
B{t)A{z)a{t -z) + {■■■), D<^^\z)B'-^^t) = a<^^\z - t)B<'^^t)D<-^^z)R'-^^\z -t) + {■■■). Between two 
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commutations of A{z) with one B{t), we push a factor a{t — z) to the right, and between two commu- 
tations of D{z) with one B{t), we push a factor R to the right, using the relations (|l7|). Once A{z) and 
D{z) are completely to the right, they act on v g{X) according to ( |l3| , [l^ ). We get eventually 



a^"^ (tm - z, Ai - 7, a) a("-i) (u-i - z, Ai - 7, A + jh^"''>) 



... a( 



/CI 

i)(ii-z,Ai-7,A + 7^A(j") 



tr 



(0) 



fc=2 

a(°)(2-^i,Ai -^7,A) a(°H2-^2,Ai - (£+1)7, A) 



7(1'-"HAi-7,A) 



. . . a(o) (z - Un , Ai - + m - 1)7, A) (z) i?^'"") (z - A 
^(o,m-i)|^^_i^_^^A + 7/i(™)) ... 7?("'i)(z-ii,A + 7^;iW 



fc=2 



(A). 



What we have in mind is to be sent back to the diagonalization of a transfer matrix T of the same 
kind as T, with dimension N decreased by 1. This is the case if we choose the right dependence of the 
vector g in the variable Ai. 

Proposition. If we take g^i- - '™)(A) G(A) ip{X) with 
(18) 

where ci is an arbitrary constant, then 
(19) 



N e 



G(A) = e^'^^Y[Yl0iX^-X,-jn), 

j=2 p=l 



gt-'^'H^) = G(A) 



S+{z;ti, . . . ,t,n) e 



9iz) 



where S±(z;ti,. .. , t^) = 



e{z - Ne-f) 
e{z-tk±j) 



S^{z;t,,... ,t,^) t(™.-.i)(z) 



and 



p(l,... ,m) j-(m,... ^-p(l,...,m)^ ^ 



tr 



(0) 



m— I 



(0) 



[/p io a conjugation by f — (f^^^ • . . f^™)) (where f is the shift operator T = d\a,g{T i)i^2,... ,n ) , 

f is the transfer matrix acting in W — (C^-i)'8(w-i)£ giygji })y expression similar to Ma Uw. 
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Proof. Let us look at the first term of go, which is given by the action of a diagonal matrix on g. 
Because 17 is a zero weight vector, wc find that this term is simply a multiple of g given by 



5(1— ™)(A-7C^l) 



N 



n u{z-tk+-f) -fj 6>(Ai-Aj--7) (1 
oi^-t.) n ,(A,-A,-(£+ih) ■ - 

For this to be exactly g'^^'- - '™)(A) up to a constant factor, we have to take 

5(1'- ^")(A) = G(A) Vi(A) 

with the function G as in (p^. 

Now let us simplify the second term. If we compute the jth element of the matrix a(z — ti, Ai 
£7, X) ■ ■ ■a{z - tm, Xi - (i + m - 1)7, X)^{z), we find 

[a{z - ti, Ai ~£j,X) ■■■ a{z ~ t^, Xi - {£ + m - 1)7, A) 

0(z) 0{z - tk - 7) 0{Xi - X, - h) ^ 



(z - (m + e)-f) 



n 



(z - tk) 9iXi - Aj) 



and then it is clear that 



[a{z - ti, Ai - £7, A) • • • a{z - Ai - (£ + m - 1)7, A) $(z)]^. G(A) 
0(z) "0(z-ffc-7) 



n 



G(A) r,. 



So when we take the trace we simply get the transfer matrix 7" up to a scalar coefficient. □ 

The first "unwanted term" is obtained similarly, except that for the commutation of A{z) and D(z) 
accross B{ti), we keep the second part of the commutation relations: 



g['---"'\X) = 



m 

i^"^(t™-ti,Ai-7,A) ••• a(2)(t2-ii,Ai-7,A + 7^/iW) 



-z,Ai -7,A + 7^/iW 



(1,... .m)(\ _ 



(Ai-7,A) 



tr 



(0) 



c(o) (z-h,Xi- £7, A) a^''' ih -t2,Xi-{£+ 1)7, A) 



,(0) 



,(0) 



(ii-<™,Ai-(£ + m-l)7,A) i?("")(<i-i„,A) •■ 



•• i?(°'^)(ii-t2,A + 7E^^'0 



fe=3 
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(l,...,m) 



(A). 



Proposition. With 5(i—™)(A) = G(A) V'(A) where G(A) is given by (Uqj, we /ia?;e 



0{h - N£j) 



eih-z) 



^(A) 
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where 



k=2 



id X is the diagonal {N — 1) x {N — 1) matrix of elements [X{z,t, X)]jj = 



e{z-t + Xj - Ai + e-f) 



^" 9{Xj - Ai + £-f) 

forj=2,...,N. 

We shall not compute directly the other "unwanted terms" . To obtain them, let us change the order 
of the B operators in (15), using the commutation relations 17). We get 

(20) ^{t, , . . . , t„. A) - 5(2) . . . (t„)i3(i) (t,) y 7^(l^™-■ ^'") (A), 

where 

rn 

^(i;™,...,2) ^ 7?(i™)(ii-t„,A) •■• i?(i2)(ti-i2,A + 75]/i('^). 

Starting with (^0|) and applying the same procedure as before, we find another expression for the 
"wanted term"; let us check that the two results are equivalent. We find now 

-1 



9^0 



S+{z;t2, . . . ,im,ii) e ''"^ 

^(l;m,...,2) 



-"'(A) = G(A) (7^(l^™^■■■,2)y 

Since S± are symmetric functions of (ii, . . . , im), and thanks to some commutation property of T and 
'R. which is written below, it is straightforward to see that this second expression is equal to ([l9|). 

Proposition. 

(21) j-(l,m,...,2)(-^-j ^(l;m,...,2) ^ -^(l;m,... ,2) ^(m,... ,1) (^^^ 

Proof. We can write T^'"— = t^o) l(0;m.,... ,i) ^^^^i 

£(0;m,... ,1) ^ (f(ni)'j j^(Om) (m-1)^ ^(Om-1) . . . ^^^^ ^(01)p (0)f (1) . . . f(m) 
^ p(0)^(Om) ^p(m)^ ^(Om-1) (m-l)^ . . . J^(01)p (2) . . . p(m) 

because f ("•'f commutes with R^^^\ Using the dynamical Yang-Baxter equation (|l|) in the form 

p(o)^(oi) ^fWY^B.'-"^^ (t^^^Y^ ^''^^^ = (f(^''^~^ ^(°^''f(°)^/"^) ' 
one can see that £(";i'2)-^(i;2) ^ ^(i:2)£(0;2,i) ^ rp^^^j^ ^^^^^ relations 

^(l;m+l:... ,2) _ ^(lm+1) ^p(m+l)^ ^(l;m,... ,2)p(m+l) 

^ (^'^0^^ |^p(m + l)^"^ ^(lm+l)p(l)^(l;m,....2)p(m+l)^ 

it is easy to prove the result by recursion on to. □ 

Let us note that at this point in the usual nested Bethe ansatz, the situation is a little simpler. Indeed 
if there is no shift entering the definition of the Lax operator, by cyclicity of the trace T'^^'™" - '^'(z) = 
T(™--'i)(z), and 7t(i;'".-- .2) jg precisely equal to T^i-™'- - '2) (ii). In that case, the relation ^) is just 
the expression of the commutation of the transfer matrices for different values of the spectral parameter. 
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But here because of the dynamical feature of the Lax operator, 7^(i'™> -- ,2) j-^^^- ^ transfer matrix, and 
^(i,m,...,2)(^-) is not equal to T^™-- ■i)(z). 

It is now clear that the expression of the second "unwanted term" is: 

,(^--)(A) = G(A) 



X(2)(z,t2,A) 



eit2 - Z) 



9{t2 - Ne-f) 

which, with the use of relation (|2l|), can be written 

0(7) 



92 



(l,...,m) 



(A) = G(A) 



0{t2-z) 



1 be written 



^(l;m,...,2) 



Oit2) 



0{t2 - Nl-i) 



V'(A). 



The expressions of the other "unwanted terms" are obtained by repeating several times the permu- 
tation of the B's. 

Proposition. The cancellation of all the "unwanted terms" is equivalent to the set of relations 

r('"'--i)(i,)^(A) = e--- Ojt.-Nij) ^ e{t,-t,+^) 



B{tk) 



Oytk -ti - J) 



6. The Bethe equations 

After the first step of the Bethe ansatz described in the preceding section, we are led to the prob- 
lem of the diagonalization of T^^^^^^^'--- '^\z) for arbitrary values of the parameters ti. Let us call 
^((^-i)^.--- ,1) the zero weight eigenstate of T^'^^"^^^'- - (z), with eigenvalue e{z) (both e and if) 
depending implicitely on the parameters ti). If the following conditions are satisfied: 

6[tk) fj^ 0(tk-ti~j) 



i=l 



the vector V'l^ii • ■ ■ i i(N-i)ti ^) given by ( |15[ ) will be an eigenstate of T{z) with eigenvalue 
, . ^ -pr 0{z-tk+l) ejz) yr e{z-t,-j) 

^ ' e{z-tk) + 0{z-m^) li e{z~tk) '^'>- 

The second step of the Bethe ansatz thus consists in diagonalizing Tii^-^)^'--- '^^(z) by repeating the 
same procedure. The total shift f (i>- - '(^-1)^) entering the definition of T*^*^^"^^^'- - '^^ has no action 
because we are looking for zero weight eigenstates only. We seek eigenstates in the form 

where v is the vector of W given by w = 62 ■ ■ • (8) 62 and Bj = L^j for j = 3, . . . , iV. Everything works 
out similarly, except the fact that the operators Dij ~ dj for i, j = 3, . . . act a little differently on 
V because of the parameters tk in T: 

N 

b{z) [g(\) v] = ^e,j $j(z,A) g(A-Lj,) w, 

i=3 



12 
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where^, z,A = ' \ J] nr , y 

We continue in this way until the A''th step; at this last step, the diagonalization problem to solve is 
simply i/)(AAr — 7) =Sn tpi^N), the solution of which is (up to a constant) i^iX^) = e'^"'^'^ , £jv = e~''"^'^. 
The whole procedure can be summarized as follows. 

Definition. We introduce a set of parameters . . . , i(jv_„)^) for n = 0, . . . ,N — 1, with 

tf = {j-ih, wj = i,...,m. 

We define a set of functions Sn{z), n = 1,. . . ,N, by the recursion relation 

{N-n)f: . („) , (N-n+l)e . ,(n-l)x (N-n)l . („) , 

where ci, . . . ,cn are arbitrary complex parameters. 

For each n = l,... ,N — 1, we define a Lax matrix by setting: 



N N N 

« \ . \ . ^. . fCA ^^^^ 



R (n\z, A) = ^ en (g) eu + ^ a{z, Xi - Xj) en ® Cjj + ^ I3{z, Aj - A^-) ® Cj 



'(.V-n + l)£ 



where the product on k is taken in increasing order if n is odd, and decreasing order if n is even. The 
notation j txi k means j > k if n is odd, and j < k if n is even. T{n) is the matrix diag(ri)i=„^... _jv 
The creation operators are the elements of C given by 

Bj{n\t) = Cnj{n\t), j = n+l,...,N, 

where the subscripts nj refer to the auxiliary space (0). 
Finally we define the states ipniXn, ■ . ■ ,Xn) by il^NiXw 

) = e'=-'>'^-^' and the recursion relation (for 

n=l,... ,N -1) 
where 



AT I: 

Gn{Xn,... ,Xn) = e^"^" H H OiXn-Xj-tn)- 

j=n+l p=l 

Proposition. If the parameters are solution of the Bethe equations 

(■Y'O^' ^(4") -4") +7) o{tt'-tt-'^-,) ^(^^-r^)) 

l\ ^(4") _ _ ^) 11 e{tt^-tt-'^) fi ^(4")_4"+i) + ^) 

37^ Is 



e 



7(cn-c„+i)^ Vn = l,... ,iV-l, Vfc= 1,... ,{N -n)i, 
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then ■0i(Ai, . . . , Ajv) is a zero weight eigenstate ofT{z) — tr'-*'-'£ with eigenvalue 

e,(z)=;^e-- n ^ r n n ^ r n ^ ^ 



'i=l )n=l 9{z~tm) \ „-^^i 9{z tm) m=l 0{z — tm '—7) 

Proof . The cancellation of all the "unwanted terms" is equivalent to the set of equations 

' ^ ~ l\ eit^^t!^-^) il ^(4") -4") -7)' 

for n = 1, . . . , iV — 1 and k — 1, . . . , (A^ — n)i. It is easy to see that 

£„(zj-2^e YV o(,_t(^)) 11 11 11 e(z^t^'-'^ 

i=n m=l '^\^ '-m ) j=n \ m=l ''l^ '-m ) m=l '^{'^ l) 

and so in the expression of £n+i(4"^) there is only one non zero term, corresponding to i = n + 1. □ 
Remark. Since the Ruijsenaars operator M is related to the transfer matrix T{z) by the relation (p^, 



"01 (Ai, ... , Aat) is an eigenfunction of M with eigenvalue e = ^1/ Z'^f) ^li^)- This quantity does not 



depend on z, so we can evaluate it at z = 0. We find 
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